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AN ALL-PURPOSE COMPUTER PROGRAM FOR LIDAR SIMULATIONS 
AND REDUCTION OF LIDAR DATA 

Charles Braun, National Environmental Satellite and Data 
Information Service, NOAA, Washington, DC, 20233, USA 1 

A general computer program has been written for the dual purpose of 
simulating lidar return signals and extracting H 2 0 vapor and temperature (T) 
profiles from actual lidar returns. This program will be used to retrieve H 2 0 
vapor and T profiles from data obtained with the University of Maryland 
ground-based lidar system. 

The program is designed to simulate and to analyze lidar returns from a 
monostatic system for two cases: when the output laser pulse is monochro- 

matic, and when the pulse has a known wavelength spectrum. For the purpose of 
simulation studies, the backscattering from the atmosphere is assumed to have 
two components, an elastic part due to backscattering from aerosols and an 
inelastic part due to Rayleigh and Brillouin scattering. The simulation part 
of the program has several model H 2 0 vapor, T, and aerosol profiles which can 
be combined to form different model atmospheres. The simulation program also 
permits the testing of different methods of numerical filtering to suppresss 
noise, and testing the accuracy of different methods of numerical integration. 

Initially the data reduction portion of the program will be checked by 
using it to retrieve H 2 0 vapor profiles from lidar signals generated in the 
simulation portion of the program. It will then be used to recover H 2 0 pro- 
files from lidar data obtained in the lower troposphere (z=0 to 5km). 
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The algorithm for 1^0 retrieval is based on an iterative method developed 
by Zuev et al. (1983) for use in the lower troposphere. In the lower tropo- 
sphere, vapor absorption linewidths are large enough compared to Doppler 
widths so that Rayleigh scattering can be treated as elastic scattering. The 
ratio of on-line to off-line signals to be analyzed is then given by 
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z = altitude 


p„ = HpO vapor density 
H^u <- 

K = HgO vapor absorption coefficient 


fi (v-v ) = normalized frequency spectrum for 

n laser pulse as a function of frequency 

v-v measured from center of 
absorption line. 

Some numerical and graphical examples of simulated and retrieved profiles 
will be given to demonstrate the operation of this program. 
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